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In  the  application  of  marine  cables  as  tension  members,  serious  structural  failures  have  occurred 
as  a result  of  hockling  (forming  of  a loop  under  torsion  combined  with  insufficient  tension).  This 
report  endeavors  to  establish  a criterion  on  end  torques  and  forces  for  preventing  hockling.  It  is 
shown  that  for  rods  in  tension  the  torque  corresponding  to  Greenhill’s  100-year-old  formula  de- 
fining the  onset  of  bending  is  also  the  maximum  torque  which  can  be  sustained  at  that  tension. 
Thus,  a straight  rod  or  cable  loaded  by  tension  and  the  Greenhill  torque  is  unstable,  and  Green- 
hill’s  formula  modified  by  an  appropriate  safety  factor  provides  the  desired  no-hockie  criterion.  . 
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20.  Abstract  (Continued) 


For  long  straight  rods  or  cables  under  tension  T and  twisting  moment  M,  this  criterion  for  stabil- 
ity is  < 4TEI,  where  El  is  the  bending  stiffness  under  load.  The  approach  used  was  to 
numerically  solve  in  nondimensional  form  the  two-point  boundary-value  problem  of  the  rod  under 
axial  end  torques  and  forces.  The  report  includes  curves  showing  torques  and  forces  for  the 
possible  range  of  deflection  curves. 
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GREENHILL’S  FORMULA  AND  THE  MECHANICS 
OF  CABLE  HOCKLING 


INTRODUCTION 

In  the  application  of  marine  cables  as  tension  members,  such  as  in  lifting  objects 
from  the  ocean  floor,  serious  structural  failures  have  occurred  as  a result  of  two  phenom- 
ena called  hockling  and  bird  caging.  Both  problems  stem  from  the  torsional  moments 
which  develop  in  the  cable,  usually  as  a result  of  its  tendency  to  unwind  under  load. 
Hockling  occurs  when,  for  a given  torque  reaction  at  the  ends  of  the  cable,  the  tension 
becomes  insufficient  to  keep  the  cable  taut,  resulting  in  the  formation  of  a loop  or 
hockle.  When  increased  tension  is  subsequently  reapplied  to  a hockled  cable,  the  loop 
tends  to  tighten,  causing  the  cable  to  fail.  Bird  caging,  on  the  other  hand,  occurs  when 
individual  wires  or  strands  unravel  under  excessive  torsional  loads  that  are  reverse  to  the 
direction  in  which  the  wire  strands  are  wound. 

The  hockling  of  cables  would  seem  to  be  closely  related  to  the  well-known  problem 
of  the  elastica  (thin  elastic  rod  under  combined  tension  and  moments  applied  at  its  ends). 
This  problem  has  received  much  attention  in  the  literature.  Kirchhoff  (1,2]  recognized 
that  the  deflection  curve  of  the  elastica  is  governed  by  the  same  set  of  differential  equa- 
tions as  the  motion  of  a heavy  spinning  top.  Greenhill  [2-4]  gave  a buckling  formula 
for  the  rod  subje  'ted  fo  tension  (or  thrust,  if  negative)  and  twisting  couples,  based  on  the 
assumption  of  in  in:.c  >s  mal  bending  deformations.  Southwell  [5] , in  a discussion  of  the 
elastica  under  ena  1c  rc  i with  zero  moments,  showed  that  a column  under  thrust  at  and 
above  the  Euler  buckling  load  remains  stable.  More  recently,  considerable  effort  has  gone 
into  studies  of  the  dynamic  characteristics  of  thin  three-dimensional  beams.  An  excellent 
review  and  bibliography  of  this  subject  through  1972  may  be  found  in  Ref.  6.  The  prob- 
lem of  rods  of  variable  cross  section  has  been  comprehensively  addressed  by  Green,  Naghdi, 
and  Wenner  [7] . 

In  spite  of  all  this  attention  to  the  problem  of  the  elastica,  some  puzzling  questions 
have  remained  with  regard  to  the  hockling  of  cables:  What  are  the  stability  characteristics 
of  small-deflection  solutions  corresponding  to  Greenhill’s  formula  with  respect  to  increasing 
torque  or  decreasing  tension?  What  in  fact  is  the  largest  torque  for  given  tension  to  which 
a cable  may  safely  be  subjected?  And,  conversely,  how  far  may  the  tension  safely  be 
lowered  for  given  end  torque?  When  this  report’s  author  first  became  interested  in  the 
hockling  problem,  he  did  a rudimentary  experiment,  trying  to  put  bending  deflections  into 
a straight  rod  under  twist  and  modest  tension.  The  rod  was  twisted  until  numerous  helical 
slip  lines  developed,  indicating  substantial  shear  yielding.  Yet  no  significant  bending  de- 
flections were  observed.  But  cables  in  tension  do  hockle.  Must  they  be  described  in  a way 
essentially  different  from  rods,  or  can  the  observed  difference  in  behavior  be  explained  en- 
tirely in  terms  of  the  much  lower  bending  stiffness  pertaining  to  cables? 
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The  purpose  of  the  analysis  described  in  this  report  was  to  obtain  some  insight  into 
these  questions  and  if  possible  to  establish  a loading  criterion  for  the  avoidance  of  cable 
hockles.  The  approach  was  to  obtain  computer  solutions  to  the  differential  equations  of 
the  elastica  in  nondimensional  form  for  the  full  range  of  values  of  axial  end  forces  and 
moments. 


STATEMENT  OF  THE  PROBLEM 

The  problem  is  to  calculate  the  deflection  curve  and  strain  energy  of  a cable  or  rod 
that  is  prismatic  and  straight  when  unloaded  and  is  “hinged”  at  each  end  so  as  to  permit 
only  an  axial  moment  M and  an  axial  force  T.  Axial  means  along  the  straight  line  joining 
the  two  end  points.  Although  this  presents  a two-point  boundary-value  problem,  it  can 
be  solved  in  reverse  as  an  initial-value  problem  by  assuming  a starting  angle  y between  the 
deflection  curve  and  the  axis  and  then  solving  the  differential  equations  step  by  step  un- 
til a point  possessing  appropriate  symmetry  properties  is  reached.  This  point  (or  alterna- 
tively the  end  point)  must  be  established  from  the  local  properties  of  the  deflection  curve, 
since  the  length  is  initially  unknown  and  must  be  calculated.  When  the  problem  is  solved 
in  nondimensional  form,  a two-parameter  set  of  solutions  is  sufficient  to  cover  the  full 
range  of  possible  loading  values. 


DERIVATION  OF  EQUATIONS 

A length  V.  of  rod  or  cable  that  is  straight  and  prismatic  when  unloaded  and  possesses 
a bending  stiffness  El  and  torsional  rigidity  GJ  is  assumed  to  lie  along  a space  curve  and 
is  referred  to  a fixed  coordinate  system  OXYZ  with  unit  vectors  I,  J,  and  K as  shown  in 
Fig.  1.  The  +X  axis  is  vertical  upward.  Arc  length  along  the  cable  is  denoted  b>  s and 
increases  as  shown  in  the  figure.  The  cable  is  assumed  to  be  “hinged”  at  O through  some- 
thing resembling  a universal  joint  such  that  the  reaction  there  consists  solely  of  a vertical 
force  T and  a vertical  couple  M,  both  directed  downward  if  positive.  T and  M denote  the 
magnitudes  of  T and  M (0)  respectively.  At  the  free  end  P of  the  cable,  where  x(s)  = 

[x,  y,  z],  the  required  equilibrating  forces  are  a constant  vertical  upward  force  T and  a 
couple  M (s),  where 


T(s)  = T ( 0)  = T = T\ 


(la) 


and 


M(s)  = M (0)  + T X x(s)  = Ml  - TzJ  + 7>K. 


(lb) 


At  the  free  end  P of  the  cable  segment  OP,  Fig.  1 shows  the  unit  vectors  of  two 
other  coordinate  systems:  £,  tj,  J and  t,  n,  b.  Here  f = t is  the  unit  tangent  vector  to 
the  space  curve  defined  by  the  cable’s  center  line.  The  vectors  £ and  tj  are  perpendicular 
to  the  tangent  and  thus  lie  in  the  plane  of  the  cross  section  of  the  cable.  They  are  taken 
as  principal  axes  of  the  cross  section  and  are  considered  to  be  inscribed  in  the  cross  sec- 
tion (rotate  with  the  cable). 


' i ih  ■|- 
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On  the  other  hand.  .-*  ««««»  «*  P™-  P™*1  “d  «“  ““ 1 °’ 
the  center  line. 


I,  i,  usual  181  to  introduce  an  “angular  velocity”  vector  O «»  denote,  the  vector 
Iate  of  t»  ol  the  system  (•  U,  t «th  respect  to  arc  length  a. 

«™  . i.  . unit  vector,  its  arclength  derivative  is  thus  given  by 


dt/ds  = t - S2  X t. 

Cross  multiplication  by  t yields 

fi  = Kb  + Sift, 


where  n and  b are  defined  by 


and 


t’  = SI  X t = Kn 


b = t X n 
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and  Q.b  = K is  the  principal  curvature  of  the  center  line.  In  terms  of  the  local  compon- 
ents of  the  moment  vector  M the  components  of  H are 


n,  = M,  /OJ, 


nn  = Mn  = 0, 


nb  = Mb/EI  = K. 

Differentiation  of  (2)  further  yields  the  derivative  of  the  normal  vector, 


Taking  components  of  (4)  in  the  t and  b directions  respectively  gives 


n'  ■ t = -K 


n • b = -b  1 n.  (5) 

The  quantities  n'  • n and  b'  • b on  the  other  hand  are  zero,  since  they  represent  the  com- 
ponent of  the  derivative  of  a unit  vector  upon  itself.  The  remaining  required  derivative 
component  is  b'  ■ t,  which  is  zero  as  shown  by  the  following:  Since  t ■ b = 0,  it  follows 
that  t'  ‘ b + t • b’  = 0.  But  since  t'  is  perpendicular  to  b,  then  t’  • b and  hence  also  t ■ b' 
must  vanish.  The  preceding  calculation  of  the  components  of  the  t',  n',  b'  vectors  along 
the  n,  b coordinates  may  be  summarized  as  the  matrix  equation  representing  the 
Frenet-Serret  differential  equations  for  a space  curve: 


0 K 0 t 


n = -K  0 if  n (6) 

b1  0 0 b 

The  quantity  ip,  substituted  for  n'  • b in  (5),  is  thus  seen  to  be  the  geometric  torsion, 
or  tortuosity  as  Love  [2]  calls  it,  of  the  space  curve  traced  by  the  center  line  of  the  cable 
or  rod.  Differentiating  equation  (lb),  noting  that  x'  = t and  using  equations  (3),  gives  along 
the  t,  n,  b components 
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M'  • t = M',  = 0, 


K(<pEl  - M, ) + T • b = 0., 


EIK'  = -T  n.  (7c) 

Equation  (7a)  states  that  the  internal  twisting  moment,  and  hence  the  twist  angle  per  unit 
length,  is  constant  along  the  arc  length.  Mt  = Mt(0)  is  one  of  the  three  first  integrals 
available  to  this  problem,  analogously  with  the  problem  of  the  heavy  top.  Equation  (7b) 
is  needed  along  with  equations  (3)  and  (6)  to  complete  the  system  of  equations. 


The  following  non-dimensional  quantities  are  introduced: 


S = Ts/M,  X = Tx/M,  Y = Ty/M,  Z = Tz/M,  L = n/M ; 


F = M2/TEl,  H = KEI/M,  -1>  = <pEl/M, 


(*)  = d/dS  = (M/T)  ('),  Mt  = M//M  s o,,  (0). 

Further  the  direction  cosines,  or  components  ay  of  the  t,  n,  b vectors  with  respect  to  the 
fixed  system  I,  J,  K,  as  well  as  a set  of  Euler  angles  A,  B,  C to  remove  unncessary  re- 
dundacies  in  the  direction  cosines,  are  defined: 


an  = cos  (t,  I), 


a12  = cos  (t,  J), 


etc.,  and 


Kyi  = 


sAsBcC  - cAsC  sAsBsC  + cAcC  sAcB 


cAsBcC  + sAsC  cAsBsC  - sAcC  cAcB 


where  cB  means  cos  B,  sC  means  sin  C,  etc. 
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In  these  terms  the  complete  system  of  differential  equations  defining  the  space  curve 
and  other  quantities  pertaining  to  the  rod  or  cable  consists  of  equation  (8)  for  the  a(y  in 
terms  of  the  Euler  angles  together  with  the  following: 


H 

- 0:31  - Za32  + ^331 

(9a) 

<l> 

= Mt  - a3l/FH , if  H f0, 

(9b) 

<l> 

= Mt/2,  if  H = 0 (by  L’Hospital’s  rule). 

(9c) 

A* 

= F(HcAsB/cB  + <!>), 

(9d) 

B* 

= - FHsA , 

(9c) 

c* 

= FHcA/cB, 

(9f) 

X* 

= «n> 

(9g) 

Y* 

= al2’ 

(9h) 

Z* 

= ft13- 

(9i) 

INITIAL  CONDITIONS 

Initial  conditions  are  chosen  with  X0  = Y0  = Z0  = 0 and  so  that  the  tangent  vector 
initially  makes  an  angle  7 (0  < 7 < 180°)  with  the  X axis.  An  appropriate  set  of  initial 
values  of  the  Euler  angles  is 

= "Z2’ 


B0=0, 


Cg  ~ 7. 

thus  making  the  curve  begin  in  the  XY  plane. 

Equations  (9)  become  singular  when  cB  = 0.  This  happens  infrequently,  and  is  en- 
tirely a difficulty  of  the  particular  Euler-angle  representation.  It  is  avoided  either  by 
choosing  a different  initial  azimuth  (say  A0  = 0,  B0  = 7,  C0  = 0)  or  by  departing  from 
the  Euler-angle  formulation  whenever  cB  is  small  and  instead  using  derivatives  of  a set  of 
three  independent  direction  cosines. 


(■pwwnuB- 
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FINDING  MIDHOCKLE 


To  determine  when  midhockle,  the  point  of  symmetry  halfway  between  the  “hinged” 
ends,  is  reached,  we  seek  an  arc  length  sm  at  which  K(sm  + s)  = K(s„  - s)  and  #(s_  + s)  = 
i p(sm  - s)  for  any  s.  This  happens  when  all  the  odd  derivatives  K> 2™  + l)  (s)  and  ^(2Af+l) 
(s)  are  equal  to  zero  at  s = sm.  It  is  easily  shown  that  this  happens  whenever  K'(s)  = 0 or, 
by  equation  (7c),  when  a2t  = 0- 

An  outline  of  the  proof  comes  in  two  parts:  The  first  part  to  be  proved  is  that  if 
K(2n  + t)  = ^(2n  + i)  = o for  n = 0,  1,  ....  Af-1  and  if  Jjf(2N+l)  = o,  then  also  ^(2A7+l ) = q 
This  follows  from  successive  double  differentiations  of  (7b)  and  use  of  (6)  and  (7c).  Each 
such  doubly  differentiated  equation  consists  of  terms  which  are  products  of  derivatives  of 
order  zero  up  to  the  order  of  differentiation.  Since  the  sum  of  orders  for  each  term  is 
odd,  each  term  must  contain  at  least  one  factor  which  is  an  odd-order  derivative  of  either 
K or  . But  by  hypothesis  all  these  factors  except  the  highest  derivative  of  >p  are  zero. 
Hence  also  the  highest  derivative  must  vanish. 


The  second  part  to  be  proved  is  that  if  2,1  + 1 ) = ^(2n  + 1)  = 0 for  n 
then  also  = q.  Since  by  (7c)  differentiated  2 n times 


0,  ...,  N-l, 


P[K(2n  + l)  = _T  . n(2 n)( 

it  suffices  to  show  that  -T  • nl2JVl=  0,  or  that  n(2JV)  has  a zero  I component.  Let  the 
direction-cosine  matrix  of  (8)  be  denoted  by  r and  the  curvature-torsion  matrix  of  (6)  by 

x: 


and 


X = 


Then  (6)  may  be  written  as 


0 K 

-K  0 

0 -i f 


T = XT. 


* 

0 


(10) 


To  be  shown  therefore  is  that  under  the  assumed  circumstances  the  { 2,1 } component 
of  matrix  r,2'V*  is  zero.  Successive  differentiation  and  substitution  of  (10)  shows  that  the 
higher  derivatives  of  r are  expressible  in  the  form 
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where  each  term  oi‘  fn  is  a product  of  from  one  to  n of  the  n matrices  x through  x*”’1  *• 
From  the  formai  process  of  differentiation  it  follows  that  these  terms  are  homogeneous 
of  degree  n in  each  term’s  sum  of  orders  plus  number  of  factors.  As  an  example  a term 
XX  Wx  has  u total  order  0 + 3 + l+  0 + 0 = 4 and  has  five  factors  and  therefore  could 
belong  only  to  f9.  Thus  for  n = 2N  the  highest  derivative  term  is  x2JV_1  multiplied  by  no 
other  factors.  Because  of  the  homogeneity  of  orders  plus  factors  of  degree  2N,  every  term 
of  /2jv  other  than  this  highest  derivative  term  must  also  possess  one  of  two  properties: 
Either  it  is  the  product  of  an  even  number  of  factors  or  it  must  contain  at  least  one  lower 
odd-order  derivative.  In  either  case,  the  matrix  term  can  be  shown  to  be  zero  at  all  “odd” 
locations  including  { 2,1 } ; hence  the  same  follows  for  their  sum  t^2N\ 

From  the  preceding  two  parts  of  the  proof  it  follows  that  a point  at  which  K'  = 0 
(or  a2i  = 0)  is  always  a point  of  symmetry. 


BRIEF  DESCRIPTION  OF  THE  FORTRAN  LISTINGS 

The  Fortran  routine  to  accomplish  the  required  calculations  consists  of  the  main  pro- 
gram Hockle  and  subroutines  Dcalc,  Sym,  and  Step.  The  listings  are  given  in  Appendix  A. 

The  main  program  accepts  sets  of  values  for  F,  y,  increment  of  nondimensional  arc 
length  AS,  and  print  interval  (to  permit  printing  less  than  every  calculation  step).  It  also 
calculates  initial  conditions  and  determines  when  midhockle  is  passed  (a21  changes  sign), 
at  which  time  it  interpolates  for  Sm  . Endhockle  is  prescribed  to  be  at  S = L = 2 Sm , and 
the  option  exists  to  keep  calculating  that  far.  The  program  Hockle  contains  necessary 
print  and  punch  instructions  and  calls  subroutines  as  needed.  An  important  output  is  sets 
of  values  of  nondimensional  moment  u = MV/EI  = FL  and  force  v = TV.2 /El  = FL2. 

The  subroutine  Dcalc  calculates  the  distance  of  closest  approach  Dmin  between  pairs 
of  symmetric  points  on  the  rod.  This  is  done  because  if  that  distance  became  zero  in  a 
real  cable,  then  the  mutual  lateral  forces  at  the  point  of  intersection  would  no  longer  per- 
mit representation  by  the  present  mathematical  model,  which  would  allow  the  cable  to 
pass  right  through  itscl.'  The  calculation  of  Umm  is  performed  using  points  only  between 
zero  and  midhockle  and  does  not  require  stepping  past  midhockle. 

The  subroutine  Sym  does  require  calculation  all  the  way  to  endhockle  and  calculates 
the  degree  of  symmetry  about  midhockle  of  selected  quantities  such  as  <l>  and  a21.  This 
provides  a check  on  the  accuracy  of  calculation,  since  in  theory  the  symmetry  should  be 
perfect. 

Finally,  the  subroutine  Step  performs  the  basic  stepping  procedure  on  the  differential 
equations  and  calculates  the  required  quantities:  direction  cosines,  curvature  H (both  di- 
rectly and  from  a first  integral  for  comparison),  >l> , etc.  When  informed  by  Hockle  that 
cB  is  too  small  to  use  the  Euler  angles,  Step  uses  alternate  equations  for  stepping.  The 
stepping  interval  is  the  inputted  AS  except  either  side  of  midhockle,  where  steps  and 
print  points  are  chosen  in  such  a way  as  to  preserve  symmetry  of  the  calculations  and 

print  intervals.  [ 
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DISCUSSION  OF  RESULTS 

Relation  Between  End  Forces,  End  Moments, 
and  Cable  Configuration 

Once  the  program  has  calculated  the  nondimensiona)  length  L of  a symmetric  space 
curve  for  given  values  of  load  parameter  F and  initial  angle  y,  then  the  values  of  the  non- 
dimensional  moment  and  force 


u = MV/E  / = FL 


and 


v = TV2 /El  = FL 2 


are  readily  calculated.  The  results  are  shown  in  Fig.  2 for  a number  of  values  of  y,  with 
a few  F = constant  curves  also  indicated.  The  y = 0°  and  y = 180°  curves  constitute  the 
two  halves  of  a parabola  corresponding  to  Greenhill’s  formula 


v = 1 (u/2)2  - it2  I , (11) 

referred  to  in  the  Introduction.  The  reason  for  the  absolute-value  signs  and  thus  the  para- 
bola with  the  lower  portion  reflected  in  the  u axis  is  that  there  is  the  following  duality  of 
solutions:  If  a (F,  y)  input  yields  a («,  e)  solution,  then  a (-F,  n-y)  input  yields  the  iden- 
tical (u,  - v ) solution.  However  the  concept  of  positive  F and  T representing  a tension  and 
negative  F and  T representing  a compression  makes  sense  only  for  the  Greenhill  case  of 
y = 0 (or,  for  reversed  signs,  180°).  For  example,  the  elastica,  under  end  force  only,  can 
go  with  a continuously  increasing  force  from  Euler  buckling  in  compression,  at  (0,  it2)  in 
the  figure,  to  a looped  rod  in  tension,  with  no  clear  break  in  between.  To  exhibit  this  be- 
havior correctly,  all  forces  as  well  as  moments  are  therefore  shown  in  the  first  quadrant, 
with  the  result  that  half  of  the  Greenhill  parabola  is  reflected  in  the  u axis.  Also,  because 
of  these  same  symmetries,  all  y = constant  curves  are  perpendicular  to  the  v axis  at  u = 0. 
In  addition  the  y = 90°  curve  is  perpendicular  to  the  u axis  at  v = 0. 

An  important  set  of  check  points  is  provided  on  the  v axis.  Clearly  this  corresponds 
to  the  elastica  with  zero  moment.  As  shown  by  Southwell  [5] , the  relation  on  u - 0 be- 
tween v and  y is  given  by 


v = 4K2(fr). 

where 


It  = sin|(7T-y)/2]  and 

where  R(k)  is  the  complete  elliptic  integral  of  the  first  kind.  This  formula  checks  the 
computer  solution  for  small  values  of  F and  hence  of  u. 
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Fifi.  2 — Nomlimi'nsional  end  moments  and  forces  on  an  elastica  for 
various  end  ancles  y and  load  parameters  F -M  ^ IT  FA  or  u~lv 


In  Fig.  2 all  nontrivial  solutions  to  this  problem  as  expressed  in  the  uv  quarter  plane 
fail  between  the  two  segments  of  the  Greenhili  parabola.  (Trivial  solutions  exist  through- 
out the  uv  quarter  plane.)  All  combinations  of  moment  and  force  inside  the  180°  seg- 
ment of  the  parabola  correspond  to  the  trivial  solution  of  twist  only,  and  zero  bending. 
Similarly  moments  outside  the  0°  segment  are  unattainable,  since  only  an  unstable  trivial 
solution  exists  there. 

The  dashed  curve  Dmm/L  = 0 along  with  the  portion  of  the  e axis  above  21.55  is 
the  locus  in  uv  space  along  which  the  calculated  Reflection  curve  of  the  loaded  cable  pos- 
sesses a self-intersection.  This  locus  provides  a barrier  to  the  configuration  of  a real  cable. 
For  example,  if  a cable  loaded  by  a constant  thrust  v = 20  is  subjected  to  an  increasing 
moment  u,  a self-intersection  will  occur  when  u equais  approximately  3,  and  any  config- 
uration corresponding  to  larger  moments  will  not  be  correctly  predicted  by  the  present 
theory.  With  reference  again  to  Southwell’s  solution,  the  Dmm/L  = 0 curve  meets  the  v 
axis  when  the  ratio  of  the  complete  elliptic  integral  of  the  second  kind  i«(h)  to  the  one  of 
the  first  kind  equals  1/2: 


£(ft)/K(fc)  = 1/2- 

This  in  turn  occurs  when  v = 21.55  and  y - 49.29°  as  shown  in  Fig.  2. 
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If  a rod  is  subjected  to  a constant  compressive  force  below  the  Euler  load  while  the 
twisting  moment  is  increased  starting  at  zero,  the  rod  remains  straight  until  the  moment 
reaches  the  lower  Greenhill  value  at  y = 180°.  This  point  is  stable  however  in  that  the 
moment  can  be  increased  beyond  this  value  at  least  until  the  dashed  curve  is  reached. 

If  a rod  or  cable  begins  in  tension,  it  remains  straight  until  the  upper  Greenhill  value 
is  reached  at  y = 0°,  but  at  this  point  no  further  increase  in  torque  is  possible  and  any 
further  twisting  deformation  must  he  accompanied  by  a decreased  reaction  torque.  Thus 
for  a rod  or  cable  in  combined  tension  and  twist,  GreenhilPs  formula  does  not  represent 
merely  the  point  of  first  departure  from  a trivial  twist-onlv  solution.  It  additionally  rep- 
resents the  maximum  torque  which  can  be  applied  at  a given  tension,  or  the  minimum 
tension  necessary  to  support  a given  torque,  and  represents  a point  of  instability  if  any 
attempts  were  made  to  increase  the  torque  relative  to  the  force.  This  instability  must  be 
particularly  violent  in  the  case  of  cables,  for  the  following  reason:  For  a rod  or  cable  in 
combined  tension  and  sub-Greenhill  torsion,  all  strain  energy  is  stored  in  twist,  a relatively 
stiff  mode.  As  the  Greenhill  curve  is  approached  and  bending  becomes  possible,  a sub- 
stantial part  of  this  twisting  strain  energy  must  be  transferred  into  bending.  But  a cable 
is  distinguished  by  its  very  low  bending  stiffness,  so  that  the  conversion  of  a given  amount 
of  strain  energy  would  require  much  larger  deflections  than  would  be  the  case  for  a rod. 

The  previous  paragraphs  also  illustrate  the  double-valuedness  of  the  graph  of  Fig.  2 
in  the  area  between  the  two  branches  of  the  Greenhill  curve.  Consider  a cable  loaded  by 
a nondimensional  moment  u = 7.2  and  a nondimensional  force  v = 52.4,  at  the  point  in 
Fig.  2 at  which  y = 15°  and  F = 1.  One  possible  configuration  of  the  cable  under  this 
loading  is  the  straight  cable  in  tension,  with  y actually  equal  to  0°  and  not  15°.  This  con- 
figuration can  result  if  the  tension  on  an  initially  unloaded  cable  is  raised  to  52.4  and  a 
subsequently  applied  moment  is  raised  to  7.2.  Since  this  is  below  the  Greenhill  torque  in 
tension,  the  cable  remains  straight.  The  other  possible  configuration  under  the  same  load- 
ing is  a hockled  cable  with  y really  equal  to  15°.  This  may  be  produced  by  “compressing” 
an  initially  straight  cable  past  Euler  buckling,  forming  a loop  which  is  tightened  as  the 
force  is  increased  to  52.4  (now  a “tension”  because  the  end  points  have  passed  each  other), 
with  a subsequent  increase  in  moment  to  7.2,  applied  in  the  direction  in  which  the  loop 
is  able  to  open  partially. 

The  quarter  plane  u > 0,  u > 0 is  thus  divided  into  three  distinct  regions:  Inside  the 
180°  branch  of  the  Greenhill  curve,  only  the  trivial  solution  exists,  and  it  is  stable.  In  the 
second  region,  between  the  180°  and  0°  branches,  the  configuration  is  double  valued  for 
any  given  loading.  Mere,  if  loading  begins  from  a straight  rod  in  tension,  then  the  trivial 
s'lution  is  stable,  but  if  the  loading  begins  from  compression,  then  the  trivial  solution  is 
unstable  and  the  nontrivial  configuration  results.  In  the  third  region,  outside  the  0°  branch 
of  the  Greenhill  curve,  only  the  unstable  trivial  solution  exists,  so  that  such  loadings  cannot 
be  sustained. 

It  is  instructive  to  examine  from  an  energy  point  of  view  the  possible  rod  or  cable 
configurations  at  the  load  point  in  Fig.  2 (u  = 2 n,  v = 0)  where  all  the  y = constant  curves 
intersect.  The  loading  corresponding  to  this  point  supports  equilibrium  configurations  for 
all  values  of  y from  0°  to  180°.  The  deflection  curve  corresponding  to  these  configura- 
tions would  in  general  be  a complete  turn  of  a helix  of  pitch  rr /2  - y.  It  is  a helix  because 
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the  force  is  zero  and  the  internal  moment  therefore  is  constant;  it  contains  a full  revolu- 
tion because  the  end  moments  act  along  the  direction  of  the  line  joining  the  end  points. 
For  y = 0°  or  180°,  the  helix  is  reduced  to  a straight  line  in  pure  twist  under  end  mo- 
ment A/;  for  y = 90°,  it  is  a full  circular  hoop  in  pure  bending  under  end  moments  M 
perpendicular  to  the  plane  of  the  circle.  Of  all  possible  equilibrium  configurations,  the 
stable  one  is  that  one  for  which  the  applied  end  moments  are  at  the  lowest  possible  level 
of  potential  energy,  i.e.,  the  one  for  which  the  strain  energy  is  maximum. 


For  this  case  of  pure  moment  loading,  the  total  strain  energy  is 


1 

U = Ml  U/2EI  + Mf  K/2 GJ, 


where  GJ  is  the  torsional  stiffness.  In  terms  of  the  applied  moment  M and  the  angle  y, 
Mb  = M sin  y and  Mt=M  cos  y , so  that 


U = JW2t/2  (sin2  y/EI  + cos2  y/GJ ). 
The  extreme  values  of  U occur  when  dU/dy  = 0,  or  when 


sin  2y  (1/EI  - 1/GJ)  = 0, 

or  y = 0°,  90°,  and  180°.  Which  extreme  corresponds  to  the  maximum  strain  energy  then 
depends  on  the  stiffnesses,  as  may  he  seen  by  putting  the  three  values  of  y back  into  U: 


U (7  = 0°,  180°)  = M2V/2GJ,  lor  pure  twist, 


and 


U ( 7 = 90°)  = M2Q./2EI,  for  pure  bending. 

Thus  for  the  rod  or  cable  loaded  by  end  moments  MV/EI  = 2rr,  the  stable  configuration 
is  the  straight  twisted  ono  if  GJ  < El  and  is  the  circular  hoop  if  GJ  > El.  For  typical 
solid  rods,  GJ  * 0.7 El,  while  for  cables  GJ  >>  El.  Thus  a hinged  cable  loaded  with 
u = 0 and  « increasing  from  zero  will  remain  straight  until  u = 2n,  when  it  will  snap  into 
the  shape  of  a closed  circular  hoop. 

In  spite  of  this  apparent  difference  between  the  configuration  of  cables  and  rods  at 
( u = 2 ir,  i = 0),  it  should  be  borne  in  mind  that  this  point  does  in  any  case  represent  the 
limit  of  stability  of  the  straight  form,  for  rods  as  well  as  for  cables.  This  behavior  of 
rods  is  also  discussed  in  Ref.  2 (paragraph  272(d),  page  417),  where  the  more  general 
Greenhill  formula  is  also  developed. 

The  main  purpose  in  presenting  Fig.  2 is  to  show  the  central  importance  of  Green- 
hill’s  formula  in  evaluating  the  stability  of  rods  and  cables  and  to  clarify  itt  relationship 


12 


NRL  REPORT  7940 


I 


F 


I 


to  the  force-only  solution  of  Southwell.  The  figure  may  also  be  used  to  determine  the 
possible  configuration  (or  configurations,  in  the  double-valued  region  between  the  two 
branches  of  the  Oreenhill  parabola  ) when  end  loads  u and  u to  the  left  of  the  7=0 
branch  are  given.  In  this  application  the  barrier  locus  Dmm/L  = 0 must  be  taken  into 
account. 


Why  Don’t  Rods  or  Solid  Wires  Hockle? 

Understanding  the  fundamental  importance  of  Greenhill’s  formula,  we  can  now  eval- 
uate under  what  conditions  a rod  in  combined  tension  and  torsion  might  undergo  bending. 
For  a circular  steel  rod  or  wire  with  Young’s  modulus  of  E = 2 X 10]1  N/m 2 and  shear 
yield  strength  os  = 350  X 106  N/m2,  the  precondition  for  bending  is  that  any  tension  be 
so  low  as  to  keep  the  actual  tensile  stress  below  0,15  X 106  N/m2,  or  a mere  0.027%  of 
a tensile  capability  of  say  550  X 106  N/m2!  Any  attempt  to  induce  bending  due  to  twist 
at  a higher  applied  tension  will  cause  torsion  shear  failure  instead.  This  conclusion  is  ar- 
rived at  as  follows:  Assume  the  tension  T is  such  that  bending  is  just  possible  by  Green- 
hill’s  formula  while  at  the  same  time  the  shear  yield  limit  us  = 2M/na2  is  reached,  where 
a is  the  radius  of  the  rod.  For  an  infinitely  long  rod  (which  would  bend  most  readily), 
Greenhill’s  formula  reduces  to  M = 2(7T7)U2  = a2  (TEn)^2.  In  the  equation  us.  = 2A//7T03 , if 
a2(TEn)112  is  substituted  for  M and  then  o,  (the  tensile  stress  produced)  is  substituted  for 
T/na2,  the  result  yields 


a,  = os2/4£. 


For  the  values  assumed  earlier  this  becomes  o,  = 0.15  X 106  N/m2  as  stated. 


EXPERIMENTAL  VERIFICATION 
Woods  Hole,  MIT,  and  NCEL  Tests 

Field  and  laboratory  tests  to  investigate  the  hockling  or  “kinking”  properties  of 
oceanographic  cable  have  been  conducted  by  Berteaux  and  Walden  [9],  Vachon  [10],  and 
Liu  [11,12],  The  laboratory  tests  by  Berteaux  and  Walden  concentrated  on  measuring 
rotation  rather  than  torque  and  therefore  cannot  be  evaluated  here  without  further  knowl- 
edge of  the  torsional  rigidity. 

Vachon  measured  for  a number  of  cables  and  loadings  the  bending  stiffness  El  (which 
he  was  also  able  to  calculate  to  good  accuracy),  the  force  7',  and  the  hockling  torque  M 
both  positive  (tightening)  and  negative  (unwinding).  His  17-foot  cables  were  long  enough 
to  reduce  end  effects  to  a modest  10%  or  so  (which  can  be  estimated  from  Greenhill’s 
formula,  equation  (11),  by  dividing  through  by  l’2  and  comparing  n2/V2  to  T/EI  or 
(M/2EI)2).  This  order  of  magnitude  is  borne  out  also  by  his  evaluation  of  the  effects  of 
end  mounting  on  hockling  torque. 

When  these  end  effects  are  ignored,  by  putting  t!  = °°  in  the  Greenhill  formula,  and 
each  experiment  is  assigned  a Greenhill  number  G = M2  /4TEI , which  theoretically  should 
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then  he  equal  to  1,  the  results  show  G to  range  from  about  0.5  to  5.  These  values  of  G 
are  at  least  the  right  order  of  magnitude  and  have  a reasonable  mean  value.  Vachon’s 
results  show  a rather  wide  divergence  in  positive  and  negative  hockling  torques  for  iden- 
tical cables  and  tensions.  If  the  torque  had  been  consistently  higher  in  the  positive  di- 
rection, one  might  be  tempted  to  ascribe  this  divergence  to  increased  bending  stiffness 
for  a tightened  cable.  Unfortunately,  in  half  the  cases  the  divergence  goes  the  other  way, 
and  no  explanation  can  be  offered  at  this  time.  In  the  one  case  in  which  the  positive  and 
negative  hockling  torques  fell  to  within  about  10%  of  their  average,  the  calculated  Green- 
hill  numbers  turn  out  to  be  an  encouraging  1.2  and  0.8  respectively. 

Reading  about  the  explosiveness  with  which  some  of  Vachon’s  kinking  experiments 
were  completed  (“a  3/16  inch  thick  circular  ring,  supporting  the  water  can  [used  as  a 
weight],  was  straightened  out,  thus  dropping  the  can  to  the  floor”)  is  more  amusing  than 
having  observed  it.  But  in  retrospect,  this  explosive  behavior  is  quite  consistent  with  the 
instability  of  the  Greenhill  loading  in  tension  combined  with  the  substantial  energy  stored 
in  the  relatively  stiff  twisting  mode  at  the  instant  when  that  instability  occurs.  Vachon 
was  well  aware  of  the  energy  consideration,  which  he  reviewed  in  his  report. 

The  applicability  of  Greenhill ’s  formula  to  cable  kinking  has  also  been  recognized  by 
Liu,  whose  recent  experiments  described  in  Ref.  11  appear  to  be  well  correlated  with  the 
formula  over  a wide  loading  range.  Reference  12  is  a more  detailed  report  of  his  results. 


NRL  Tests 

In  view  of  the  incomplete  knowledge  at  NRL  about  the  conditions  under  which 
hockling  tests  at  other  Laboratories  were  conducted  and,  even  more,  to  get  some  engi- 
neering insight  into  the  phenomena  under  study,  a few  rudimentary  measurements  were  ' 
made  at  NRL.  The  most  sophisticated  measuring  tools  were  fish  scales  and  weights 
(steel,  to  keep  the  floor  dry),  and  the  time  allowed  for  the  experiment  was  about  2 days. 
Nonetheless  the  results  showed  good  agreement  with  theory. 

The  no-load  small-deflection  bending  stiffness  El  for  a sample  of  1/4-inch  cable  was 
measured  by  loading  several  lengths  of  it  as  cantilever  beams.  Calculated  values  of  El 
ranged  from  0.051  to  0.063  Nm 2,  with  0.057  as  an  average  value. 

A 12-foot  length  of  (his  cable  was  suspended  from  the  ceiling.  Because  of  a lack  of 
safety  devices,  loading  was  kept  at  light  tensions  up  to  18  N and  torques  up  to  2 Nm. 

For  four  loading  cases  Greenhill  numbers  G = M2/4TEI  of  0.97,  0.72,  0.97,  and  1.10 
were  obtained.  In  view  the  closeness  of  these  observations  to  the  theoretical  G - 1,  more 
comprehensive  and  better  instrumented  hockling  experiments  are  anticipated,  with  a view 
to  developing  a usable  method  for  specifying  torsional  properties  of  marine  cables. 


CONCLUSIONS 

The  conclusions  are  as  follows: 

• The  importance  of  Greenhill’s  100-year-old  formula  for  determining  the  elastic 
stability  of  rods  in  combined  axial  force  and  twist  far  surpasses  merely  defining  the  onset 
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of  possible  bending  modes.  For  rods  or  cables  in  tension  the  Greenhill  condition  repre- 
sents the  largest  torque  which  can  be  applied  for  a given  tension  or  the  lowest  tension 
capable  of  supporting  a given  torque.  It  represents  a point  of  instability  with  respect  to 
increasing  torque  or  decreasing  tension.  This  instability  is  expected  to  be  particularly 
violent  in  the  case  of  cables,  because  of  their  low  bending  stiffness. 

• For  applications  in  which  a cable  remains  substantially  straight,  such  as  the  lifting 
of  objects  from  the  ocean  floor,  the  Greenhill  formula,  modified  by  an  appropriate  safety 
factor,  should  provide  a valid  criterion  for  estimating  the  onset  of  hockling. 

• For  applications  involving  initially  curved  cables,  such  as  in  towing,  where  the  in- 
fluences of  gravity  and  drag  are  strong,  the  situation  is  less  clear.  It  is  known  that  the 
Greenhill  condition  cannot  be  exceeded,  but  it  is  possible  that  instabilities  might  already 
occur  at  much  lower  values  of  torque  or  higher  values  of  tension,  since  the  relevance  of 
Greenhill’s  formula  depends  on  the  relatively  large  twisting  stiffness  of  the  initially  un- 
bent and  straight  cable.  The  stability  of  initially  curved  cables  thus  is  an  interesting  sub- 
ject for  further  research. 
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Appendix  A 

FORTRAN  LISTINGS  FOR  PROGRAM  HOCKLE 


PROGRAM  HOCKLE 

DIMENSION  FF(100) ,GftM|100) , DSD (100) 

COMMON  /Sl.OKl/  A11.A12.A13.F.H.PHI.A.B.C.X.Y.Z.S.A21.A31.A32.A33. 

1 OA.DB.DC.DX.DY.DZ.aBAR.ICOUNT.HNEW.aIIO.GaMMa.IBFLaG 
COMMON  /BL0K2/  KOUN1 » SYMpHI  ( 1500) * SYMALF ( 1500 ) 

COMMON  /BL0K3/  XS(lt>00)  »'?^(1500>  #2S(XD00>  .LS.DMI N , DS* SMI N, NSTEPS 
REAL  L 

THE  FOLLOWING  FUNCTION  CALCULATES  GREENHILL  LENGTH 
XL (X) ■ 6.2831B53072  / SQRTfX  * t*-4)) 

WE  NOW  READ  IN  NUMBER  OF  STEPS  BETWEEN  PRINT  LINES  AND  STOP  SIGNAL! 
DS»0* 

001  READ  2.NSTEPS.NST0P 
2 FORMAT (215) 

IF(NSTOP.GT.O)  GO  TO  1002 

WE  NOW  READ  IN  UP  TO  100  POINlS  F, GAMMA, DS 

READ  l.NF.NG.NDS. I FF ( I ) , I = 1 , NF)  . (GAM(J) .J*1.NG> . (DSD (KOI ,KQ»1.NDS> 
1 F0RMAU3I3,  (/8F10.0) ) 

FOR  OUTPUT  PURPOSES  WE  NOW  COMPUTE  THE  LENGTH  OF  THE  PRINT  INTERVAL 
AND  INTIAUZE  THE'  CASE  nUMDER;  PRN.NCASE 

NCASEb0 

WE  RUN  THE  CALCULATION  THROUGH  THE  F, GAMMA  PAIRS, 

I IS  A COUNTER  To  STEP  I MROUSH  THE  VALUES  OF  F 
J IS  A COUNTER  TO  STEP  THROUGH  THE  VALUES  OF  GAMMA 
KQ  IS  A COUNTER  TO  STEP  THE  VALUES  OF  DS 

DO  1000  1=1, NF 
F*FF (I) 

•••  CHECK  TO  SEE  IFl  F IS  IN  THE!  INTERVAL  (0,*)  AND  IGNORE  GREENHILL! 
LENGTH  IFl  IT  IS, 

OL»0 

I r ( (F.LT.o*)  .or.  (F-GT.A.)  ) GL»XL(F) 

• **  INNER  LOOPS  ON  J ANU  Y.Q 

DO  1000  J;-i,NG 
00  1000  K C ■ • 1 , NOS 
DS=0S0 (k0! 

PRN*NSTEf:3'-DS 

IBFLA5«0 

SMnO, 

SMIN«o. 

R'O* 

RN“0 « 

NPRlNToo 

NTEST=o 

Sl«S2*S3"S4'0. 
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ABAR-0 

KOuNT“l 

MlDFLAG»0 

ICOUNT«0 

GAMMA»GAM(J) 

C 

C •••  PRINT  INITIAL!  COND11  IONS  AND  CASE  NUMBER 
C 

NCASE-NCA$E*1 
PRINT  3*NCAS£' 

3 FORMaTUHI, — INITIAL'  CONDITIONS  FOR  CASE  NO. •»!*//) 

PRINT  AiF, GAMMA 

♦ FORMAT  ( 1H  '.10Xt«F"»*F6.2.lDX*GAMMA«»,F6.2//) 

C 

r.  •*•  convert  gamma  to  rauians  Eor  calculation 

c 

GAMMA-GAMMA  ••  3 . 1 A 1 S926536/1 80  . 

C 

C INITIALIZE  angles  AND  COORDINATES 
C 

A»l. 5707963268 

B»0» 

C-GAMMA 

c 

C ALTERNATE  INITIAL  VALUES  A»o..  B-GAMMA,  C"0.  INSERT  HERE 
C 

X-0. 

r-0. 

z-o« 

S*0. 

rnmax«0. 

RN«o« 

LS«i 

c 

C CALCULATE  INITIAL!  VALUES  OF  AII.J) 

C 

AU-COS(B).COS(C)  S AllO-Ali 
ABAR-ABAR.All 
I COUNT" I COUNT* I 
A12«COS(B)»SIN(CI 
Al3=-SIN(B) 

A21"SIN(A)»sIN(B)*COsIC)  • COst A)#SlN(C) 
A3l*C0S(A)»SIN(B|»CUS(C)  ♦ SIN(A) *SIN(CI 
A3a»cOS(A)*SIN(i!)«SlN(C)  * S IN C A> «COS (C» 

A33"C0S  t A) *COS IB) 

C 

C ALTERNATE  A21  INSERT  heme 
C 

C CALCULATE  INITIAL'  VALUES  OF  H.PHI , PSI ,HNEW 
C 

H«A3W«A32*T*A33 

PhI°0'5*Aji 

IF(a0S(H!.. 000000011  5»5»Al 

*1  PHI«  All  • A31  / (F*H) 

5 hneW"Sort((2./f)»(cos(gamma)-aui  ♦ si n (gamma j »»2 > 
psi-o. 
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ITESTIbO 
IT£ST2»0 
PUNCH  60»NCASE 

60  formati*  initial  values  for  case#*I6> 

PUNCH  61,S,pHl,X,YiZ.RN»HipsI 
PUNCH  61,Ail,A12,Al3,A21,A3i,A32,A33 

61  FORMAT (8F10.6) 

66  F0RMAT1aF14.8»2FI2*B1 

c 

C print  OUTPUT  HEADERS* 

PRINT  6 

6 FORMAT (/AH  S,6X*3HPHI ,5** 1HX* 6X. 1HY» 6X, 1HZ.5X.2HH0.5X, 3HPSI »5X, 
aHA.6X.lHB,6X*lHCt6X.lHII.9X.3HAll,AX.3HA12,4X,3HA13,4X.3HA21.AXt 
23HA31 1 Ax*3HA32« AXi  3HA33/) 

IFrL.l 

PRINT  INITIAL'  VALUES 

T PRINT  8, S.PHI ,X,Y,Z,HNEN,PSI, A, B,C,H, All, A12,Al 3 , A21 , A3l , A32, A33 
B FORMATllH  ( 1 B F 7 . A ) 

RN=SQRT(Y»*2»Z»»2) 

IF(RN.GT.RNMAX)  RNmaX-RN 
NPRINTkNPRINT*! 
call  SYM(IFFL'*S1*S2*S3*SA) 

IF(MIOFLAG.EQ.I)  so  to  50 
XS(LS)«X 
YS(LS)oY 
ZS(LS).Z 
l.S-LS*i 
50  CONTINUE 

WE  NOW  BEGIN  TO  STEP  OUK  WAY  ALONG  THE  CABLE 
ISP  IS  A COUNTER  CONTROLLING  PRINT 
WE  ARE  LOOKING  FOR  MIOHUCKLE 

ISPnO 

9  A2iOLO»A2j 
CALL  STEP(OS) 

•••  TEST  FOR  MIOHOCKLE 

IF( ( (A21.A210L0) .LT.O)  .AND.  (MIDFLAG. EO. 0) > GO  TO  13 

• WE  ARE  NOT  AT  MIOHOCKLE. 

HAVE  WE  COME  TO  THE  END  OF  THE  CABLE. 

10  IF  (MIDFLAG. NE'.l)  GO  TO  11 
IFINPRINT.EQ.NTESTl  GO  TU  25 

WE  ARE  NOT  AT  END  OR  MIUHOCKLt.  CHECK  TO  SEE  IF  IT  IS  TIME  TO  PRINT 

11  13P«ISP*1 

C WE  ARE  NOW  READY  TO  PRINT  SO  WEI  NEED  TO  COMPUTE  HNEW.PSI.PHI 
PSI bASIN t Z/S0RT(Z»“2*Y**t)  ) 
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IF(ISP.LT.NSTEPS)  GO  TO  V 
GO  TO  7 

WE  ARE  AT  MIOHOCKLE 


13  PRINT  14 

14  FORMAT I////1H  ,*■ 


MIOHOCKLE 


•*///) 


WE  now  FIND  SM  The  LEN3|H  TO  MIOHOCKLE 
MIDFLAG*  1 . 

OSM«-A2i«DS/IA2l»A2lOLD) 

CALL  STEP (OSM) 

RN«S0RT(T  »‘2*Z  **2> 

IF(rN.GT*RNMAX)  rNMAX-RN 

XS(LS)«X 

YS(LS)bY 

ZS(LS)»Z 

CURRENT  VALUES  OF'  VARTAHLES  in  COMMON  are  midhockle  values  SMBS. 
PSIbaSINI  Z/SQRT (Z**2  *Y**2)  ) 

L*2.*s 

FLbF»L 

FL2»FL*L 


CALCULATE  TWIST  »NO  BENUING  ENERGIES 

ENTLbF»FL2*0.5*C0S (GAMMA) b»2 

abar«abar/icount  ........uu.ih,  ti,iiRin 

ENB»FL2"C0S  (GAMMA)  *MFL2*0 *5*SIN  ( GAMMA)  **2-FL2  ABAR 

SMbS' 

GpGAMMA«1 BO. F3. 1*15*26536 
PRINT  HEADERS  FOR  MIOHOCKLE  PRINTOUT 

17  FORmItcIH  ,32X  *PRInT*.3X.*GREENHIU,*  , 2X»* MIDHOCKLE*, 3 X,* TOTAL** 

I2*x. ‘Twist  benoing*) 

IB  FORMATC1H  ,5X,.F.,6A,*GAMHA*,4X,*STEP  SIZE  INTERVAL!  *.3  (•  LENGTH' 
1 •I.4X*FL*t&X,. FL2*,7X, -ENERGY  ENERGY*)  * 

PRINT  19.F.G  *Ds>  PRN,GL*sH, Ll*FL,FL2»ENTL*ENB 

19  FORMAT (1H*F10*5*F7,2»3X*VF10*5/F//) 


S»**F12.5l 


CALL  DCALC 
PRINT  51  OMIN*SMIN 

51  FORMAT ( 1 H *•  MINIMUM  DISTANCE  IS  *,F12.6** 

PRINT  B.S.pHl.X,Y.Z.HNEw.PSl»A»B»c*H»All*A12,Al3,A21*A31,A32,A33 

PUNCH  62 

62  FORMAT  <«HIOHOCKLE*>  „ . D1JU  » 

PUNCH  61 ,F,GAMMA,OS,PRN,GL,L,RN,RNMaX 
PUNCH  64»FL,FL2,ENTL,ENB iUMIN, SMIN 

PUNCH  61,S,PHI.*iY,Z.RNiH.PSI 
PUNCH  61 *A11*A12,A1J,A21,A31, A32*A33 

WE  NOW  GO  TO  POINT  SYMMETRIC  WITH  PRINT  VALUE  PREC.EEOING  MIOHOCKLE 
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C 


2* 


DS2e0s*DSM 
CALL  STEP (0S2> 


K2»isp 

00  24  K»1«K2 
CALL  STEP(OS) 

call  SYH(IFPU'*SliS2*S3tSA»  .... 

PRINT  B.S,PHliX,Y.ZiHNEW«PSl.A.B.C.H.AHtAl2.A13.A2l,A31,A32,A33 


CALCULATE  TO  END  OFi  CABLE  WITH  ORIGINAL  Ds 


IFFL-2 

iSPno 

SAVE  NUMBER  OF  STEPS  TO  MID  HOCKLE  FOR  END  TEST 

NTEST*NPRINT 
NPRINT»l 
GO  TO  9 

XE  GET  TO  HERE  WHEN  WE  HAVE  COME  TO  THE  END  OF  THE  CABLE 

25  IFFL'3 

CALL  SYM  < IFFLi  St iS2>S3  »S4) 

26  SYMMETRY  FACTORS*— * »6X*  *pHI».  F12.5, 1 OX  • *A21  •» 

1F1Z.5///) 

27  FORMAT  I//1h"  «•— RMS  VALUES— 15X tF12 .5*  13X*F12.5) 

R1-S1/S3 

RZ“S£''S*  . 

PRINT  28  Pi  R2 

28  FORMAT  (//]H*  t •— RATIOS**— *t  I8*»F12.5*13X»Fl2.5) 

PUNCH  63 

63  FORMAT («ENDHOCKLE*l 

PUNCH  6liS|PHI*X»YtAtRNiH*PSI 
PUNCH  61*a11*A1z«A13iAH1i«31*a32*A33 
PUNCH  61»sl»52«s3»sA*Rl’R* 

C IF  THERE  ARE  MORE1  PAIRS  p,GAMMA  TO  BE  DONE  RESET  AND  GO  AGAIN 
C OTHERWISE  STOP. 

C 

1000  CONTINUE 
GO  TO  1001 
1002  STOP 
END 
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cXoNT/0LoS3/LXSUi»OO).r5(15OO>.ZS(15OO),LS,DMIN,DS.SHIN,NSTEPS 

DIMENSION  0(15001 


THIS  SUBROUTINE  CALCULATES  RELATIVE  MINIMUM  OF  DISTANCE 
BETWEEN  PAIRS  OF  SYMMETRIC  points 


DO  5 Jal • 1500 
5 DC J» • 

0MIN»-2. 

XMaXS (LS) 

YMaYS(LS) 

ZMnZS(LS) 

MIN*LS»1 

MM»»1 

DO  10  I«liMIN 
J-LS-I 

T ■YM*#2*ZM»*2 

IF(T1 »LT. 0.0000001)  00  TO  10 
T2*(YM#ZS( J) •ZM*IS(U)  i,#2 
T3»(XM-XSf J))«»2 
D(J)*-2.»SQRTU2/Tl»l3) 

10  CONTINUE  , , . . 

print  loo*  (O(J)i  U*1»LS> 

100  FORMAT ( 1H  .10F10.5) 

c 

c LOOK  FOR  MINIMA 

c 

TEST  B0 (LS) 

I«LS"1 

z0  IF(DII).LT.TEST)  GO  TO  30 
TEST“D  ( I ) 

21  I**I”i 

IF(I.EO.O)  GO  TO  50 
GO  TO  20 

30  IF(D(l)aLT.0«l  GO  TU  21 

31  TEST =D  ( I ) 

I«Jar 

32  IF (Dll) •GT«TESTI  GO  TO  40 
TEST "0(1) 

1*1-1 
GO  TO  3? 

40  DHINbTEST 

SMIN«(  I»j)«OS»NSTEMS 
RETURN 

51  FORMATUH**  d has  no  relative,  minimum 
RETURN 
END 
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SUBROUTINE  SYM(I,S1.S2.S3»S*I 

common  /0LOKI/  *11  »M2i*13»F«HtPHI,A,B*CiX»r*Z*S*A2liA3l.A32.*33, 
1 DaiDb«DC,DX.DY.OZ.ABARi IC0UNTiHNE*,a110iGaMMa»IBFLaO 
COMMON  /DL0K2/  KOUNt , SYMPHI ( 1500 ) t SYMALF < 1 500 ) 

C 

c this  subroutine  is  called  at  each  print  step  and  calculates  symmetry 
c factors  at  end  or  cable. 
c 

GO  TO  ( 1 «2«3) •! 

C 

C #E  ARE  STILL  IN  FIRST  half  of  CABLEw  KE  ACCUMULATE  values  OF  PHI.A21 

c 

1 SYMPHI (KOUNT) "PHI 
SYMALF <KOUNT)"A21 
S3«S3*PHl*»2 
S*»S**A21»«2 
KOUNT^KOUnt*! 

J'KOUNT 

1FIKOUNT.EQ.1500)  OU  TO  5 

RETURN 

C 

C HE  ARE  IN  SECOND  half  of  cable. 

c 

2 KOUNT«KOUNT-l 

SYMPhI (KOUNT)«(SYMPnI (K0UNT)-PHI)»«2 
SYMALF ( KOUnT) ■ (SYMaLF (KOUNT ) *A21 ) **2 
S1«S1»SYMPHI(K0UNT) 

S2«s2*SYMALF(K0uNT) 

RETURN 


NE  ARE  AT  THE  END  OF  THt  CABLE, 


Sl«SQRT<Sl  /(J*ll 
S2*SQRT ($2/1 J*l) ) 
S3«SQRT(S3/(J«1)I 
Sa»SQRT(Sa/(JMI) 

RETURN 
PRINT  A 
FORMAT  (lH 
STOP 
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SUBROUTINE  STEP (DS> 

COMMON  /8U0KI/  All»*12«A)3.F*H»PHIf A,B.C.X,r,Z,Si*21#*31.A32,*33, 

1 DAtDBiDC,DX»DriOZ.ABAR«ICOUNT*HNEW,AUO.GAMMA»I«FUQ 
COMMON  /BL0K2/  KOUNT i SVMHMI ( 1500) tSYMALF ( 1500) 
c this  SUBROUTINE  STEPS  The  CALCULATION  from  POINT  to  POINT  AL0N3  the 
C CABLE. 

C INPUT  DS  ■ STEP  SIZE  FOR  CALCULATION 

C 

c 

C OUTPUT  A * 


C B -ANGLES 

C c • 

C X ♦ 

C Y ‘COORDINATES 

C z ♦ 

c 


C CHECK  IBFLAG, 

C I8FLAGM  IF  B TOO  CLOSE  TO  90  DEGREES 

C IF  B TOO  CLOSE  TO  9u  WE  USE  ALTERNATE  FORMULAE 
C 

IF(IBFLAG.EQ.l)  go  10  S00 

c 

c-»«  first  we  calculate  increments  for  angles 
c 

400  OB*  (-F  • H • SlN(AI)  • Ds 

DA*  F • (H  • COS (A)  • SIN1B) /COS (B)  ♦ PHI)  • DS 
DC*  ( F * H • COSUI/COSIB))  »OS 
C 

C •••  WE  ADD  INCREMENTS  TU  PRESENT  VALUES  TO  OBTAIN  VALUES 
C AT  A POINT  DS-UNITS  FURTHER  ALONG  THE  CABLE. 

C 

A-A.DA 

B=B«D8 

C»C*DC 

SlGNCB*SlGN ( 1 « » COS (B) ) 

C 

C CHECK  To  SEE  IF  B TOO  CLOSE  TO  90  DEGrEEs 
C 

CB*ABS (COS (B) 1 
IF (CB.LT »0.1  ) OU  TO  10Q 

C 

C •*  WE  ARE  NOW  REAOY  TO  CALCULATE  INCREMENTS  TO  COORDINATES 
C 

DX-AU  • Os 
DY.A12  • OS 
0Z»At3  * OS 
C 

C •••  WE  AOO  THE  INCREMENTS  TO  THE  PREVIOUS  VALUES  OF  THE  COORDINATES 
C TO  OBTAIN  VALUES  AT  A POINT  OS-UNITS  FURTHER  ALONG  THE  CABLE. 

C 

x*x*ox 

Y*Y‘DY 

Z*Z‘QZ 

c 

C FINALLY. WE  CALCULATE  NEW  CABLE  LENGTH 
C 
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S«S*05 

•••  CALCULATE1  NEW  VALUED  FOR  A21 #A31 *A32*A33 

A2l»SlN(A)«SIN<B)*CUSIC>  - COS ( A) «SIN C C> 
A31-C0S(A>*SIN(B)»COS(C>  ♦ SIN<A)#SIN(C) 
A32«C0S(A)»SIN(B)»S1N(C)  » SlN ( A) «C0S ( C> 
A33»C0S(A)»C0S(B) 

•«*  calculate  nek  values  for  mi»ai2>ai3  and  return 


AU«C0S(B)«C0S(C) 

A12.C0S(B)«S1N(C) 

A13=»S!N(B) 

300  H«  AAl-Z»A32*Y»A33 

abarSabar*aii  s icount=icount*i 

PHI«0.5»A110 

IMABS(H).. 00000001)  16,16,15 
15  PHI»A110-A3l/(F»H) 
lfc  HNEW-0. 

HNEw2«l2./F) •< COS (GAMMA) -All) * S1 N( GAMMA) **2 
IF(HNEW2.QT.0.)HNE*»SaRT(HNEW2) 

RETURN 

GET  TO  HERE  WHEN  B TOO  CLOSE  FOR  FIRST  TIME. 

WE  SAVE  CURRENT  VALUES  OF  AiB.C 
WE  THEN  SET  IBFLA3.lt  A=C«0. 


100  AGAVE«A 
BSAVEsB 
CSAVE“C 

IBFLAGbj 

A»o. 

C«0. 

B«o* 

NOW  CALCULATE  A22tA23  FOR  FIRST  TIME  THROUGH 

A22«sIN(ASAvE) .SIN(BSAVE)*SIN(CSAVE).C0S(ASAVE)*C0S(CSAVE9 
A23nSlN(ASAVEJ«C0S(bSAVE) 

C 

C WE  NOW  CALCULATE  HEW  VALUES  FUR  XtTiZtS 

c 

200  F»X*A11»DS 
Y«Y*a12«0S 
Z«Z*A13*0s 

s.s.os 

c 

C NOW  CALCULATE'  NEW  VALUES  OF  AdtJ) 

C 

DAllrF»H<‘Acl*Os 

DA12»F6H»A?2*0S 

0*21  «F*  < — M «•  A 1 1 *PHl»A3J  )#DS 

DAl3StF«H<>A23»0S 

DASlS'-F^PHI^AZl^OS 

0A32S=-F<'PHl"A22»DS 
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I Jlf  *1'  Ill'll  .11  M 


FELIX  ROSENTHAL 


A11«A11»DA|1 

A12*A12.DA12 

A2l«A21*0A21 

Ai3S»Ai3*0Ai3S 

Al3*S0ST(i.-Ali»»2-»l2»«2) 

IF  l(A13*A13S)  .L’T.O. ) Al3«*A13 
A3 ls« A31 *DA3 1 S 

A3lBS0RT(l.-AU.»2.A21»»2) 

IF(tA3l*A3lS>.LT.O.)  A31"“ *31 

A23«  ( «A  1 1 *A  1 3#Aiil*Al2»A3l)/tl»*All*#2) 

A22«(A12»A23-A3i)/Al3 

A32S«A32*0a32S 

A32=SQRT(i.-A12«»2*A22»«2) 

IF  C ( A32*A3pS» .LT.O* » A32*"A32 
A33«(An»A31*A12»A32J/t-A13) 

GO  TO  300 


C— — — — 

c 

C GET  TO  HERE  IF  IBF|.AG«1 
C CHECK  TO  SEE  IF  B STILL1  TOO  CLOSE. 

C 

$00  CB«SQRT(1.-A13*»2> 

IFICB.LT. 6.1  ) 80  TO  2?0 

GET  TO  HERE  IF  B NO  LONGER  TOO  CLOSE* 
*E  RESTORE  B TO  PROPER  UUADRANT 

IFISIGNC8.GT.0.I  C8«-CB 
SB--A13 

DETERMINE'  BiA.C  UP  TO  ROTATIONS 

B*ACOS(CB) 

IF (SB.LT  0.)  Bb-B 
CAb  A33/CB 
SA«A23/CB 
A°ACOS (CA) 

IFISA.LT.O.)  A»-A 
CC«All/C8 
SC-A12/CB 
C“AC0S(CCI 

ifisc.lt. 6.)  C«-C 

c 

c NOW  DETERMINE'  ROTATIONS 

c 

BR0T«(B-BSAVE5/6.2B3lB530T 
NbBROTaO  »5 
B=B.n»6, 283105307 
CR0Tb(C-CSAVE1  /6. 283165307 
N=CROT*o«5 

C»C«N*6.283l85307 
AROTb(A»ASaVE'J/6.ZBJ1B5307 
N»ArOT»0.5 
A«A*Ne&. 283185307 
C 

c CONTINUE  WITH  ORIGINAL  FORMULAE! 


IBFLAG«0 
GO  TO  A00 
END 
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